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PREFACE 
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1.0  INTRODUCTION 

The  turbulent  mixing  of  two  streams  inside  a duct  occurs  in  many  devices  of 
interest  to  the  engineer.  Typical  examples  are  jet  pumps,  composite  propulsion  systems 
such  as  air-augmented  rockets,  and  various  low-speed  combustors.  The  flow  field  within 
engine  test  cells  often  involves  the  mixing  of  a high-speed  exhaust  jet  with  a low-speed 
secondary  flow.  In  many  applications,  one  must  consider  the  occurrence  of  exothermic 
chemical  reactions  in  the  mixing  process.  And  some  flows,  with  low-speed  secondary 
streams,  have  regions  of  recirculating  flow  along  the  duct  wall. 

Ducted  flows  without  recirculation  (Fig.  la)  can  be  accurately  predicted  by  solving 
the  boundary-layer  equations,  in  which  lateral  pressure  gradients  are  neglected  and  the 
velocity  gradients  in  the  lateral  direction  are  assumed  to  be  much  larger  than  those  in  the 
axial  direction.  After  the  specification  of  a semi-empirical  model  for  the  turbulent 
transport  of  momentum,  energy,  and  mass,  the  boundary-layer  equations  can  be  solved 
by  finite-difference  methods  (e.g..  Ref.  1)  or  by  integral  methods  (e.g.,  Refs.  2 - 6). 

In  the  duct  flow  field  without  recirculation,  the  axial  velocity  components  are 
directed  downstream  throughout  the  entire  flow.  However,  if  the  initial  velocity  of  the 
secondary  stream  is  very  low  compared  with  that  of  the  primary  stream,  then  a 
recirculation  zone  can  occur  along  the  duct  wall  (Fig.  lb).  If  the  entire  secondary  flow 
is  entrained  by  the  mixing  layer  before  the  mixing  layer  propagates  to  the  duct  wall,  then 
recirculation  must  occur  to  satisfy  the  requirement  for  additional  mixing  layer 
entrainment. 

Experimental  investigations,  such  as  that  of  BarchUon  and  Curtet  (Ref.  7),  have 
identified  the  main  features  of  ducted  flows  with  recirculation.  The  wall  static  pressure 
gradually  rises  between  the  duct  entrance  and  the  front  stagnation  point  (FSP),  which  is 
the  station  at  which  recirculation  begins  (Fig.  lb).  The  wall  static  pressure  then  remains 
nearly  constant  between  the  FSP  and  the  station  where  the  outer  edge  of  the  shear  layer 
reaches  the  wall  (Note  that  the  edges  of  the  shear  layer  are  somewhat  arbitrarily  defined). 
Downstream  of  this  station,  which  is  the  beginning  of  the  reattachment  region,  the  wall 
static  pressure  rises  rapidly.  In  addition,  the  reverse-flow  velocities  along  the  wall,  which 
are  largest  at  the  onset  of  reattachment,  decrease  quickly  in  the  reattachment  region  and 
become  zero  at  the  rear  stagnation  point  (RSP).  Downstream  of  the  RSP,  the  axial 
velocity  components  are  everywhere  positive,  and  the  flow  eventually  becomes  a fully 
developed  pipe  flow. 

The  experiments  indicate  that  the  reversed  flow,  between  the  edge  of  the  shear  layer 
and  the  duct  wall,  is  highly  turbulent  and  well  stirred.  Therefore,  the  flow  in  this  region 
is  nearly  one-dimensional,  with  small  gradients  in  the  axial  velocity  and  species  profiles. 
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In  essence,  the  flow  between  the  axial  stations  defined  by  the  FSP  and  the  onset  of 
reattachment  appears  to  be  a free  turbulent  jet  that  is  mixing  with  a nearly 
one-dimensional  outer  flow.  Of  course,  the  velocity  of  this  outer  flow  varies  considerably 
in  the  axial  direction. 


The  wall  boundary  layers  are  very  thin  in  the  region  of  recirculating  flow,  and  wall 
viscous  effects  seem  to  be  important  only  far  downstream  of  the  RSP  (Ref.  4). 

Flows  with  recirculation  are  truly  elliptic  in  nature  because,  in  certain  regions,  the 
axial  flow  gradients  become  large  compared  with  the  lateral  gradients.  Therefore,  an 
accurate  analysis  of  recirculating  flows  must  be  based  on  solution  of  the  Navier-Stokes 
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equations.  During  the  past  several  years,  much  progress  has  been  made  on  numerical 
techniques  for  solving  these  equations.  For  example,  the  finite-difference  analysis  of 
Gosman,  et  al.  (Ref.  8),  has  been  widely  used  for  predicting  recirculating  flow  fields.  But 
in  spite  of  the  potential  of  the  finite-difference  approach,  there  are  several  reasons  why 
this  approach  cannot  yet  be  considered  a reliable  analytical  tool  for  routine  engineering 
applications  (Refs.  9 and  10).  First,  with  an  adequately  fine  mesh  for  the  flow 
configurations  of  interest,  the  computation  time  is  large.  Second,  there  are  unresolved 
accuracy  problems  in  the  numerical  algorithms;  finite-difference  solutions  of  combustor 
flows  commonly  do  not  conserve  fuel  species  to  an  acceptable  degree.  Third,  specification 
of  boundary  conditions  is  difficult  and  is  often  done  in  an  arbitrary  manner.  Finally,  as 
with  any  analysis  of  turbulent  flow,  an  adequate  semi-empirical  model  of  the  turbulent 
transport  processes  must  be  specified;  with  the  finite-difference  approach,  the  transport 
model  must  be  applicable  to  every  point  in  the  flow  field. 

There  is  a simpler  approach  to  the  engineering  computation  of  ducted  flows  with 
recirculation.  The  integral  form  of  the  boundary-layer  equations  can  be  used  to  describe 
flows  with  recirculation  if  the  profile  shapes  for  velocity,  species,  etc.,  are  specified.  Hill 
(Refs.  4 and  11)  used  this  integral  approach  to  predict  constant-density  ducted  flows 
with  recirculation.  Of  course,  the  streamwise  turbulent  transport  terms  are  neglected  in 
this  approach,  which  results  in  a poor  description  of  the  flow  in  some  regions.  But  the 
overall  flow  field  predictions  are  adequate  for  many  engineering  purposes,  and  the 
integral  method  is  computationally  efficient.  In  addition,  the  very  nature  of  the  integral 
method  ensures  that  momentum,  species,  etc.,  are  conserved. 

In  this  report,  the  integral  approach  is  applied  to  the  prediction  of  variable-density 
ducted  turbulent  flows  with  and  without  recirculation.  Flows  with  equilibrium  chemical 
reactions  are  included  in  the  formulation.  The  flow  geometry  considered  is  either 
axisymmetric  or  planar,  and  provision  is  made  to  calculate  flows  with  an  arbitrary 
centerbody.  In  the  "basic"  integral  analysis  described  in  Section  2.0,  the  streamwise  • 
turbulent  transport  terms  are  neglected.  The  predictions  of  the  basic  analysis  are 
compared  with  available  experimental  results  in  Section  3.0,  and  the  regions  of  the  flow 
that  are  poorly  predicted  are  identified.  An  "extended"  integral  analysis,  in  which  the 
streamwise  turbulent  transport  terms  are  included,  is  described  in  Section  4.0. 

2.0  FORMULATION  OF  BASIC  INTEGRAL  ANALYSIS 

The  mathematical  approach  used  in  the  present  analysis  is  similar  to  that  used  for 
the  "l-D  Core  Theory"  of  Ref.  2.  Indeed,  the  present  analysis  can  be  considered  an 
extension  and  a refinement  of  the  "l-D  Core  Theory,"  which  is  applicable  only  to 
axisymmetric  flows  without  recirculation. 
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2.1  FLOW  REGIMES 

Several  flow  regimes*  can  occur  in  a ducted  turbulent  flow.  These  regimes  are 
characterized  by  the  width  and  position  of  the  turbulent  shear  layer  and  by  whether  or 
not  the  flow  recirculates.  The  integral  analysis  must  be  formulated  differently  for  each  of 
the  flow  regimes.  The  flow  regimes  that  can  occur  in  a nonrecirculating  flow  are 
illustrated  in  Fig.  2,  and  the  regimes  that  can  occur  in  a ducted  flow  with  recirculation 
are  illustrated  in  Fig.  3.  Thus,  there  arc  six  possible  flow  regimes  and  five  possible 
sequences  of  these  regimes. 


a.  Mixing  zone  propagates  to  centerbody  before  it  propagates 
to  duct  wall 


b.  Mixing  zone  propagates  to  duct  wall  before  it  propagates 
to  centerbody 

Figure  2.  Possible  sequences  of  regimes  for  flow 
without  recirculation. 


Mn  the  nunicnclature  used  to  describe  the  regimes,  "first"  indicates  a regime  in  which  there  is  a 
potential  cure  of  jet  fluid  between  the  centerbody  and  the  inner  edge  of  the  shear  layer.  "Second” 
indicates  a regime  in  which  the  shear  layer  extends  to  the  centerbody  but  has  not  yet  reached  the  duct 
wall.  "Wake"  indicates  a regime  in  which  a region  of  reversed  flow  exists  between  the  outer  edge  of  the 
shear  layer  and  the  duct  wall. 
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Note  that,  when  the  initial  value  of  the  secondary  velocity  is  zero,  the  "first"  regime 
(Fig.  3)  will  not  occur,  and  recirculation  will  begin  at  the  duct  entrance  (in  the  "first 
wake"  regime). 


b.  Recirculation  begins  in  first  regime  and  mixing  zone  propagates 
to  duct  wall  before  it  propagates  to  centerbody 


c.  Recirculation  begins  in  first  regime  and  mixing  zone  propagates 
to  centerbody  before  it  propagates  to  duct  wall 
Figure  3.  Possible  sequences  of  regimes  for  flow  with  recirculation. 


2.2  BASIC  ASSUMPTIONS 

The  principal  assumptions  used  in  the  analysis  are: 

1.  The  flow  is  steady  and  either  planar  or  axisymmetric. 
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2.  All  gases  obey  the  perfect  gas  law. 

3.  The  flow  is  described  by  the  integral  form  of  the  boundary-layer  equations. 

4.  The  turbulent  Prandtl  and  Schmidt  numbers  are  unity. 

5.  The  mixing  layer  is  fully  turbulent  (negligible  laminar  transport),  and  the  initial 
boundary  layers  are  negligible. 

6.  The  mixing  layer  velocity  profiles  are  shape-similar  and  are  represented  by  a 
cosine  function. 

7.  The  inviscid  portions  of  the  primary  and  secondary  flows  are  one-dimensional 
and  isentropic  (first,  first  wake,  first  wall,  and  second  regimes). 

8.  For  mixing  with  simultaneous  chemical  reactions,  the  reactions  are  in 
equilibrium. 

9.  The  duct  and  centerbody  surfaces  are  frictionless  and  adiabatic. 

10.  The  static  pressure  is  constant  in  the  first  wake  and  second  wake  regimes. 

11.  In  the  wake  region,  between  the  outer  edge  of  the  shear  layer  and  the  duct 
wall,  the  axial  component  of  velocity  and  the  element  species  are 
one-dimensional  (first  wake  and  second  wake  regimes). 

2.3  BASIC  EQUATIONS 

The  basic  boundary-layer  equations  are: 

Continuity 


^ {pur“)  + (pvr“)  = 0 
5*  dr 


where 

a = 0 for  planar  flow 
a = 1 for  axisymmctric  flow 

Axial  Momentum 


pur®  ^ F pvr®  ( pfr®  — r®  iLli 

dr  \ dr  / 5* 

where  e is  the  turbulent  eddy  viscosity. 


(1) 


(2) 
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Energy 


r7x 


pvr«|Il  . ^ (p.r“|li) 
or  dr  \ or  f 


(3) 


Element  Species  Conservation 

piir®  di;  \ p\x^  dl-  = {p<x°-  dL^ 

(9\  f9r  f7r  V dx  ) 


(4) 


Nomenclature  for  the  integral  analysis  is  shown  in  Figs.  4 and  5.  Equations  (1),  (2), 
and  (4)  are  integrated  (Ref.  2)  to  obtain  five  integral  equations:  (1)  a continuity  equation 
for  the  entire  flow,  (2)  a momentum  equation  for  the  entire  flow,  (3)  a momentum 
equation  for  the  flow  between  the  centerbody  surface  (re)  and  rm  (where  = ri+  b/2),  (4) 
a species  conservation  equation  for  the  entire  flow,  and  (5)  a species  conservation  equation 
for  the  flow  between  re  and  r^ . 


a.  Mixing  zone  propagates  to  centerbody  before  it  propagates 
to  duct  wall 


b.  Mixing  zone  propagates  to  duct  wall  before  it  propagates 
to  centerbody 

Figure  4.  Nomenclature  for  integral  analysis  of  flow 
without  recirculation. 
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Figure  5.  Nomenclature  for  integral  analysis'  of  flow 
with  recirculation. 

The  integral  equation  are: 

Overall  Continuity 


w 

/ ^ W '■“‘I''  = Pc“c‘‘o 


dr 


— our  ZL 

= dx  "■  dx 


(5) 


C 

where  the  subscript  "c"  indicates  properties  along  the  centerbody  surface  and  the  subscript 
"w"  indicates  properties  along  the  duct  wall. 


Overall  Momentum 


f (pu^)r“dr  = p « n ^ r“  '*'^w  _ [?  w ~ c 1 “P 

/ <?*  ^ ^ dx  dx  (^TT) d7 


(6) 


Half-Radius  Momentum 


f (pu^)r“dr  — u f _d  {pQ)r°dr  = r r“  - (u  ■-  u ) p u r“  - l-*^  m ~ c 1 

J dx  ^ J dx  " dx  (a  - 1)  dx 


(7) 


where  the  subscript  "m"  indicates  the  properties  along  rm  and  Tm  is  a turbulent  shear 
stress. 


Overall  Species 


/w 

-d.  (puC)  r®  dr  = p u C r®  ‘‘c  — p u C r® 
P Pc  c c c Pw  «■  w w 


(8) 


Half-Radius  Species 


mm  , 

f ^ (paC)r®dr  - J ^ (pu)r®dr  = q^r®  - p,u^(C^  - C„,)r«.^  (9) 
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where  qm  is  a turbulent  species  flux. 


Because  the  flow  is  one-dimensional  for  rc  < r < n and  for  (n  + b)  < r < rw , the 
integral  equations  can  be  rewritten  as 

Overall  Continuity 

r.— b 

ra+l  pO-HT  frO+l  fr 

- -J-  ^ J / ^ ^ (pw“w^ 

(a  -f-  1 ) n X a X (a  + I ) ax 


dr 


= o u r** £ — n u r“  y 

Pc  C C Pw  W W J.. 


dr 

n 

dx 


(10) 


r.+  b 


Overall  Momentum 

- (■'  _i  (pa^jr-dr  + i_  (p , ) 

(all)  » (3x  (a  + 1)  dx 

I 

u^r“  ^ - p r®  ^ 


,a+i  _ ^a+l 

W C 

(a  + 1) 


1 ^ 

dx 


Half-Radius  Momentum 


r,a+i  a+n  „ 

J»J G — JjL  (p_u^)  I f (pu^)r®dr  - u J=_! £ — J _d_  (p„u  ) 

(a  1-1)  dx  “(9x  (a  1-1)  dx 

i 

/*"  Ar  rr0t+^  - rO+n  J 

(pu)r®dr  = r r®  + (u  — u ) p u r“  — £ — L -£ J _E. 


r — b 


(11) 


(12) 


Overall  Species 

r a-H  ,a-t-n  /•'  , fr®"'''  - (r  -t- b)®~n  , 


p u C r®  £ — p u C r® 

Pc  C C C Pw  W W W (J3J 


Half-Radius  Species 

rr®+*_r®^n  r'”  , fr^+l  _ r®+n 

-Lj G — i jL  (p,.u  C j -I-  f (puC)r®dr  - C li-J £ — i (p  u 

(a  -I-  1)  dx  ' ^ dx  (a  -I-  1)  dx 

i 

r 

m I 

- C f (pu)r“dr  = q r“  - p u (C  — C ) r“  G. 

mjj  'r'  >min  Pcc  c " 


dr, 

dx  ' ■™  "■  ■ ^ ^ dx'  (14) 


Equations  (10)  through  (14)  apply  to  all  six  regimes  shown  in  Figs.  2 and  3.  Of 
course,  r,  goes  to  in  the  second,  second  wake,  and  third  regimes.  In  the  first,  first 
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wake,  and  first  wall  regimes,  Q = Cj  = 1,  and  Pc  and  u^;  are  the  isentropic  variables  (pj  and 
Uj).  In  the  first  and  second  regimes,  = 0,  and  Pw  and  are  the  isentropic  variables  (pa 
and  Ua).  The  pressure  gradient  terms  are  zero  in  the  first  wake  and  second  wake  regimes. 


2.4  SHEAR  LAYER  PROFILES 

The  shear  layer  velocity  profile  is  described  by 


U — U I , 

!!L  = I - L cos 

u — 11  2 2 


(15) 


The  velocity  at  the  half-radius  control  surface  (rm ) is 


= 4 “ J 


(16) 


For  a constant-pressure  nonrecirculating  flow,  the  element  species  concentration 
profile  is  related  to  the  velocity  profile  by 


C — C u — u 

1 


c.-c, 


U — U 
J ^ 


(17) 


Similarly,  the  total  enthalpy  profile  is  related  to  the  velocity  profile  by 


H - H 

w 

H . - H 

j ^ 


11  . — u 
J ^ 


(18) 


Thus,  C,  H,  and  u are  linearly  related  in  a constant-pressure  nonrecirculating  flow.  This  is 
the  well-known  Crocco  integral  solution  obtained  from  Eqs.  (2),  (3),  and  (4).  In  the 
ducted  mixing  problem,  however,  the  pressure  gradient  terms  are  negligible  only  in  the 
first  wake  and  second  wake  regimes.  In  addition,  u^,  C^,  and  H„  vary  with  x,  even  in 
the  constant-pressure  wake  regimes.  Therefore,  Eqs.  (17)  and  (18)  cannot  be  used.  The 
approach  that  has  been  taken  in  this  study  is  to  assume  that  the  concentration-velocity 
relation  can  be  expressed  as 


(19) 

where  K can  be  a function  of  x.  The  mixing  zone  concentration  profiles  described  by 
Eqs.  (15)  and  (19)  are  shown  in  Fig.  6 for  various  values  of  K.  Note  that  a power-law 
relation  similar  to  Eq.  (19)  has  been  used  to  describe  profiles  that  deviate  from  the 
Crocco  relation  because  of  nonunity  transport  coefficient  ratios  in  free  turbulent  flows 
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(Ref.  12).  Also  note  that  the  linear  relationship  between  H and  C is  retained,  even  for 
flows  with  recirculation. 


2.5  COMPUTATION  OF  DENSITY 

For  unity  Prandtl  and  Schmidt  numbers,  the  stagnation  enthalpy  and  the  gas 
properties  are  uniquely  related  to  C,  and  the  density  field  is  calculated  with  the 
techniques  described  in  Ref.  2.  Namely,  the  stagnation  temperature,  specific  heat,  and  gas 
constant  are  specified  as  a function  of  C,  which  varies  from  zero  (pure  outer  stream  gas) 
to  one  (pure  jet  gas).  For  chemically  frozen  flows,  only  the  To,  Cp,  and  R values  at  C = 
0 and  C = 1 need  be  specified.  For  flow  with  chemical  reactions,  an  equilibrium 
chemistry  analysis  is  used  to  specify  Tq,  Cp,  and  R as  tabular  functions  of  C. 

At  any  point  in  the  shear  layer,  u is  given  by  Eq.  (15),  C is  given  by  Eq.  (19),  and 
To,  Cp,  and  R are  determined.  The  density  can  then  be  obtained  from  the  energy 
equation  and  the  perfect  gas  law.  Therefore,  the  density  can  be  expressed  functionally  as 

p = p(u,  C,  p)  (20) 

2.6  TURBULENT  TRANSPORT  TERMS 

The  turbulent  shear  stress  (xm)  that  appears  in  Eq.  (12)  is  given  by 
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r 


m 


p f ^ 
rm  n 


(21) 


where  e is  the  turbulent  eddy  viscosity.  The  turbulent  species  flux  (q^  ) that  appears  in 
Eq.  (14)  is  given  by 


‘im 


f ^ 
’ <5r 


(22) 


The  eddy  viscosity  model  used  in  this  study  is  the  same  as  that  described  in  Ref.  2, 
Namely,  the  Prandtl  model  for  free-turbulent  flows  is  used,  with  the  coefficient  (k)  taken 
to  be  a function  of  the  half-radius  Mach  number  (Mm ). 


Thus, 


€ 


kb  U 

c 


(23) 


and 


k = [0.66  + 0.34  exp  (-3.42  (24) 

In  the  first  and  first  wake  regimes,  k,,  is  set  equal  to  0.007.  In  the  other  four  regimes,  ko 
is  set  equal  to  0.01 1. 

For  nonrecirculating  flows,  this  eddy  viscosity  model  has  been  shown  to  provide 
reasonably  accurate  predictions  as  long  as  the  secondary  velocity  (u^  i ) is  not  larger  than 
about  0.3  times  the  jet  velocity  (uj,)  (Ref  2). 

For  flows  with  recirculation,  the  eddy  viscosity  model  is  modified  somewhat.  The 
predictions  of  the  reattachment  region  (in  the  third  or  first  wall  regimes)  are  improved  if 
the  eddy  viscosity  is  held  constant,  or  "frozen,"  at  the  value  that  is  computed  at  the 
beginning  of  reattachment.  This  modification  is  discussed  further  in  Section  3.1. 

2.7  METHOD  OF  SOLUTION 

Sufficient  information  is  available  to  transform  Eqs.  (10)  through  (14)  into  a system 
of  ordinary  differential  equations  that  is  linear  in  the  derivatives  of  the  dependent 
variables.  The  transformation  procedure  is  described  in  Refs.  2 and  13. 

First  Regime  - The  flow  is  completely  described  by  four  parameters  (p,  r\,  b,  K),  so 
Eq.  (14)  is  not  used.  The  resulting  system  of  equations  is 
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-^1 


If 

dx 


N,^  + N,  + N 

^ dx  ^ dx 


dx 


N5 


(25) 


where  N = F for  the  continuity  equation,  N = G for  the  overall  momentum  equation,  N 
= H for  the  half-radius  momentum  equation,  and  N = I for  the  overall  species  equation. 
The  equations  of  the  coefficients  for  all  regimes  are  presented  in  Appendix  A.  The 
numerous  auxiliary  equations  necessary  for  computation  of  the  coefficients  are  developed 
in  the  same  manner  as  described  in  Ref.  2. 


After  the  coefficients  are  evaluated  numerically,  Eq.  (25)  is  solved  for  the  derivatives 
(dp/dx,  drj/dx,  db/dx,  dK/dx)  by  use  of  a matrix  factorization  technique.  The  resulting 
equations  for  the  derivatives  are  numerically  integrated  with  a modified  Euler  technique 
(variable  step  size). 

Second  Regime  - As  in  the  first  regime,  the  flow  is  fully  described  by  four 
parameters;  these  are  chosen  to  be  (p.  Uc,  b,  K)  and  the  resulting  system  of  equations  is: 

N,  ^ + N2  ^ + N3  ^ - N4  ^ = N5  (26) 

' dx  dx  ’ dx  ^ 


Third  Regime  - The  flow  is  fully  described  by  five  variables  (u^,  Uc,  p,  Cw,  K),  and 
Eq.  (14)  is  required.  The  resulting  system  of  equations  is 


M du 

N,  S! 

'■  dx 


Nr 


du 

c 

dx 


dx 


dC 

d X 


+ X-  djC  = N 
dx  ® 


(27) 


where  N = J for  the  half-radius  species  equation. 


First  Wake  Regime  - The  five  dependent  variables  are  (u^ , r„  b,  , K),  and  the 
resulting  system  of  equations  is 


du  .T  dr.  jL  dC 

!!  + N,  db  _ N I 

dx  dx  dx  dx 


N.  dK  = N 
^ dx  ® 


(28) 


Second  Wake  Regime  - The  five  dependent  variables  are  (u^,  u^,  b,  C^,  K),  and  the 
resulting  system  of  equations  is 


..  du 
— !!  + 


N2  ^ + N,  ^ + N.  ^ 

dx  ^ dx  ^ dx 


N,  dK  = V 
° dx  ° 


(29) 


First  Wall  Regime  - The  five  dependent  variables  are  (u*,  n,  p,  Cw,  K),  and  the 
resulting  system  of  equations  is 
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N,  ^ + N,  ^ ^ No  Ip  + N, 

^ dx  dx  ^ dx  ^ 


dC 

J 

dx 


N.  slK  = N' 

^ dx 


(30) 


Notes  on  the  Computer  Program  - Which  of  the  six  regimes  will  occur  in  a given 
flow,  and  the  sequence  of  these  regimes,  is  generally  not  known  beforehand.  Therefore,  a 
series  of  tests  is  incorporated  in  the  computer  program  to  automatically  provide  the 
correct  sequence  of  regimes. 

In  the  present  computer  program,  the  duct  wall  and  centerbody  geometries  are  each 
described  by  a pair  of  functions:  a fourth-order  polynomial  in  x describes  the  radius  up 
to  a prescribed  axial  station,  after  which  the  geometry  is  conical. 

The  basic  integral  analysis  has  been  programmed  for  numerical  solution  with  an  IBM 
370/165  digital  computer.  A typical  flow  field  computation  requires  a CPU  time  of  about 
30  sec. 


3.0  RESULTS  AND  DISCUSSION  OF  BASIC  INTEGRAL  ANALYSIS 

For  ducted  flows  without  recirculation,  the  present  analysis  yields  predictions  that 
are  nearly  identical  to  those  of  the  "1-D  Core  Theory"  of  Ref.  2.  Therefore,  the  major 
emphasis  in  this  section  will  be  placed  on  comparing  the  predictions  with  experimental 
results  for  flows  with  imbedded  recirculation  regions. 

Similarity  parameters  have  been  proposed  to  correlate  the  occurrence  and  extent  of 
recirculation  in  constant-density  axisymmetric  flow  (Refs.  14,  15,  and  16).  These 
similarity  parameters  are  applicable  to  constant-area  ducts  without  centerbodies.  Perhaps 
the  most  widely  used  similarity  parameter  is  that  proposed  by  Becker,  et  al.  (Ref.  16), 
which  is  commonly  called  the  Curtet  number: 


where  rn  is  the  primary  nozzle  radius  and 

"k  = <“jl  - “al^  (y  + “al  (31a) 

The  relationship  among  Ct,  Uai/uji,  and  is  shown  in  Fig.  7. 

Correlations  of  experimental  data  have  shown  that  recirculation  will  occur  when  Ct 
is  less  than  0.7  to  0.9.  In  addition,  experimental  data  on  the  axial  location  and  extent  of 
the  recirculation  region  in  various  duct  geometries  are  well  correlated  when  plotted  versus 
Ct. 
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Note  that  the  Curtet  number  is  not  applicable  to  variable-density  flows  or  to  flows 
in  ducts  with  variable  cross-sectional  area.  Schulz  (Ref.  9 and  17)  derived  a 
variable-density  similarity  parameter  but  found  that  it  does  not  adequately  correlate  the 
location  and  extent  of  recirculation  observed  in  variable-density  experiments.  That  is,  the 
variable-density  results  do  not  collapse  onto  the  constant-density  results  when  plotted 
versus  the  similarity  parameter. 


Figure  7.  Curtet  number  as  a function  of  velocity  ratio. 


3.1  EXPERIMENTS  OF  BARCHILON  AND  CURTET 

Barchilon  and  Curtet  (Ref.  7)  conducted  an  extensive  series  of  experiments  on 
constant-density  flows  in  a cylindrical  duct.  The  nozzle-duct  radius  ratio  (rn/rw)  was 
0.075,  and  the  secondary-to-primary  velocity  ratio  (Ugi/uji)  was  varied  from  zero  to  a 
value  high  enough  to  prevent  recirculation.  Predicted  and  experimental  results  for  all  six 
of  the  Barchilon  and  Curtet  flows  are  shown  in  Figs.  8 through  10.  Only  Test  No.  1 (Ct 
= 0.976)  was  conducted  with  a secondary  velocity  hi^  enough  to  avoid  recirculation 
along  the  duct  wall. 

Unless  otherwise  indicated  on  the  flgures,  all  the  computations  were  made  with  the 
eddy  viscosity  "frozen"  in  the  reattachment  region  (Section  2.6).  Two  of  the  flows,  Tests 
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No.  4 and  6,  were  computed  with  both  the  frozen  eddy  viscosity  model  and  the  local 
eddy  viscosity  given  by  Eq.  (23).  The  predicted  distributions  of  wall  velocity  (Figs.  8d 
and  f)  show  that  the  position  of  the  RSP  is  predicted  much  better  with  the  frozen 
viscosity.  The  need  to  arbitrarily  freeze  the  eddy  viscosity  in  the  reattachment  region  is  a 


b.  Test  Number  2,  C*  = 0.714,  Ua^/Uji  = 0.0465 
Figure  8.  Wall  velocity  distributions  for  Barchilon  and 
Curtet  experiments. 


result  of  the  very  strong  history  effects  on  the  turbulence  when  the  mean  flow  changes 
rapidly  with  distance.  Simple  local  eddy  viscosity  models  are  clearly  inadequate  for 
rapidly  changing  flows  such  as  occur  in  the  reattachment  region. 
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The  predicted  distributions  of  wall  velocity  (Figs.  8a  through  0 illustrate  a 
characteristic  feature  of  the  integral  solution.  The  magnitude  of  the  peak  reverse  flow 
velocity  is  overpredicted,  and  a cusp  is  formed  in  the  velocity  distribution  at  the 
beginning  of  reattachment.  In  the  region  near  the  mixing  zone  attachment  point, 


d.  Test  Number  4,  C,  = 0.305,  Uai/up  = 0.0170 
Figure  8.  Continued. 


streamwise  turbulent  transport  of  momentum  is  very  significant.  Of  course,  this 
streamwise  momentum  transport,  which  is  neglected  in  the  basic  integral  analysis,  tends 
to  round  the  cusp  in  the  wall  velocity  distribution. 
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In  spite  of  the  locally  unrealistic  behavior  of  the  predicted  wall  velocity  distribution 
near  the  onset  of  reattachment,  the  solution  seems  to  recover,  and  the  slope  of  the 
velocity  distribution  at  the  RSP  is  fairly  well  predicted. 


f.  Test  Number  6,  Ct  = 0.075,  Ugi/uji  = 0 
Figure  8.  Concluded. 
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Axial  distributions  of  centerline  velocity  for  the  six  Barchilon  and  Curtet  flows  are 
shown  in  Figs.  9a  through  f.  In  general,  the  predicted  distributions  are  satisfactory. 


b.  Test  Number  2,  Ct  = 0.714,  Uai/up  = 0.0465 
Figure  9.  Centerline  velocity  distributions  for  Barchilon 
and  Curtet  experiments. 
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Figures  9d  and  f illustrate  that  the  frozen  eddy  viscosity  model  yields  slightly  better 
predictions  of  the  centerline  velocity  decay  rate  than  does  the  local  viscosity  model. 


d.  Test  Number  4,  Ct  = 0.305,  Uai/uji  = 0.0170 
Figure  9.  Continued. 
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Axial  Distance,  x/r^^ 

f.  Test  Number  6,  Ct  = 0.075,  Uai/Uji  = 0 
Figure  9.  Concluded. 


Axial  distributions  of  wall  static  pressure  for  the  six  flows  are  shown  in  Figs.  10a 
through  f.  In  general,  the  predicted  wall  pressure  in  the  reattachment  region  tends  to  rise 
too  abruptly,  and  the  peak  wall  pressure  is  overpredicted;  this  overprediction  becomes 
larger  as  Ct  is  increased. 
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The  static  pressure  in  the  first  wake  and  second  wake  regimes  is  assumed  to  be 
constant.  However,  as  shown  in  Figs.  10c  through  f,  the  experimental  wall  static  pressures 


a.  Test  Number  1,  Ct  = 0.976,  Uai/uji  = 0.0648 


b.  Test  Number  2,  Ct  = 0.714,  Uai/Ujt  = 0.0465 
Figure  10.  Wall  pressure  distributions  for  Barchilon 
and  Curtet  experiments. 


decrease  slightly  with  distance  downstream  of  the  FSP,  until  a minimum  is  reached  at  a 
station  upstream  of  the  predicted  shear  layer  attachment  point. 
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d.  Test  Number  4,  = 0.305,  Uai/uj^  ? 0.0170 

Figure  10.  Continued. 
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e.  Test  Number  5,  Ct  = 0.152,  Uai/uji  = 0.0057 


f.  Test  Number  6,  Ct  = 0.075,  Ugi/uji  = 0 
Figure  10.  Concluded. 
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In  the  integral  analysis,  the  wall  frictional  effects  are  assumed  to  be  negligible.  To 
evaluate  the  validity  of  this  assumption,  a wall  shear  stress  term  was  added  to  the  right 
hand-side  of  Eq.  (11);  the  shear  stress  is  equal  to  1/2  CfwPwUw^  where  Cf^  is  the  skin 
friction  coefficient.  As  shown  in  Fig.  10a  for  Q = 0.976,  the  effect  of  even  a fairly  large 
value  of  Cfw  on  the  pressure  distribution  is  small.  And  the  effect  on  the  Uw  and  Uc 
distributions  is  negligible  (Figs.  8a  and  9a).  For  flows  with  recirculation,  the  effect  of 
wall  friction  is  even  smaller  than  that  shown  for  Q = 0.976. 

Distances  from  the  entrance  of  the  mixing  duct  to  the  front  and  rear  stagnation 
points  are  shown  in  Fig.  1 1 for  the  Barchilon  and  Curtet  flows.  The  analysis  predicts  no 
recirculation  for  Q greater  than  0.84;  on  the  other  hand,  the  experiments  show  that 
recirculation  persists  for  Q values  up  to  about  0.95.  The  predicted  effect  of  the 
secondary-primary  velocity  ratio  on  the  location  of  the  FSP  is  in  excellent  agreement 
with  the  experimental  results.  These  experiments  show  that  the  location  of  the  RSP  is 
independent  of  Q ; the  analysis  predicts  that  the  distance  to  the  RSP  is  a mild  function 
of  Q. 


Figure  11.  Extent  of  recirculation  zone  for  Barchilon  and 
Curtet  experiments. 
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The  analysis  yields  complete  spatial  distributions  of  the  axial  velocity  component  (u) 
and  the  density  (p)  in  the  flow  field.  By  integrating  the  continuity  equation  between  the 
axis  and  a particular  streamline,  the  geometry  of  that  streamline  can  be  determined. 
These  computations  were  carried  out  for  several  streamlines  in  the  recirculation  region  of 
the  Barchilon  and  Curtet  Test  No.  6 flow.  As  shown  in  Fig.  12,  the  predicted  streamline 
shapes  are  realistic;  in  addition,  they  are  in  fairly  good  agreement  with  the  experimental 
shapes  (Ref.  7).  The  largest  differences  between  the  predicted  and  experimental 
streamline  shapes  occur  near  the  duct  entrance,  for  x/ty,  less  than  one;  these  differences 
are  attributed  to  deviations  from  one-dimensionality  in  the  experimental  wake  region  u 
profiles. 


x/rw 

0 1. 0 2. 0 3. 0 


Figure  12.  Predicted  streamline  pattern  for  Barchilon  and  Curtet 
Test  No.  6 (Ua^/uji  = 0). 


Predicted  and  experimental  loci  of  u = 0 are  also  shown  in  Fig.  12.  The  predicted 
locus  is  satisfactory  for  x/r^  greater  than  one,  but  is  less  so  near  the  duct  entrance. 

3.2  EXPERIMENT  OF  HESKESTAD 

In  the  Barchilon  and  Curtet  recirculating  flow  experiments,  rn/r*  was  only  0.075, 
and  the  turbulent  mixing  layer  reached  the  axis  well  upstream  of  its  intersection  with  the 
duct  wall.  Therefore,  the  flow  pattern  was  as  shown  in  Fig.  3a  or  Fig.  3c,  and  the  entire 
reattachment  process  occurred  in  the  third  regime.  For  much  larger  values  of  rn/r^,  say 
0.5,  the  shear  layer  will  reach  the  duct  wall  before  it  reaches  the  axis,  and  the  flow 
pattern  shown  in  Fig.  3b  will  be  observed. 

Heskestad  (Ref.  18)  measured  the  wall  pressure  distribution  in  a constant-density 
axisymmetric  sudden-expansion  flow;  Uai/ujj  was  zero  and  rn/rw  was  0.498.  The 
predicted  pressure  distribution,  obtained  with  the  eddy  viscosity  frozen  in  the 
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reattachment  region,  is  in  good  agreement  with  the  experimental  distribution  (Fig.  13). 
As  shown  in  the  figure,  reattachment  begins  in  the  first  wall  regime;  however,  the 
predicted  RSP  occurs  just  downstream  of  the  transition  to  the  third  regime. 


Figure  13.  Wall  pressure  distribution  for  the  experiment  of  Heskestad. 


The  predicted  results  obtained  for  the  Heskestad  geometry,  as  well  as  those  obtained 
for  the  Barchilon  and  Curtet  geometry,  indicate  that  the  concept  of  freezing  the  eddy 
viscosity  in  the  reattachment  region  is  successful  for  a variety  of  configurations. 
Therefore,  all  of  the  calculations  described  in  the  remainder  of  this  report  were  obtained 
with  the  frozen  viscosity  model. 

3.3  EXPERIMENTS  OF  BECKER,  ET  AL. 

Becker,  Hottel,  and  Williams  (Refs.  16  and  19)  conducted  a series  of  experiments  on 
constant-density  axisymmetric  flow  with  recirculation.  The  nozzle-duct  radius  ratio 
(rn/fw)  was  0.032.  Therefore,  reattachment  occurred  in  the  third  regime.  In  addition  to 
measuring  the  mean  velocity  field  in  their  experiments,  Becker,  et  al.,  measured  optically 
the  concentration  field  by  use  of  a smoke  tracer  in  the  jet. 

The  axial  location  and  extent  of  recirculation  measured  by  Becker,  et  al.,  is  shown 
in  Fig.  14,  along  with  the  predictions  of  the  integral  theory.  Also  shown  in  the  figure  are 
the  predictions  for  the  geometry  of  Barchilon  and  Curtet;  the  predictions  for  the  two 
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geometries  differ  little  when  plotted  as  a function  of  Ct-  In  contrast  to  the  Barchilon  and 
Curtet  results  on  the  location  of  the  RSP  (Fig.  II),  the  results  of  Becker,  et  al., 
substantiate  the  predicted  trend  that  the  location  of  the  RSP  depends  somewhat  on  the 
Curtet  number. 


Figure  14.  Axial  extent  of  recirculation  for  the  experiments  of 
Becker,  et  al. 

Axial  distributions  of  jet  species  concentration  at  the  duct  wall  (C^)  are  shown  in 
Fig.  15  for  the  three  values  of  Ct  investigated  by  Becker,  et  al.  The  predicted  locations  of 
the  FSP's  are  indicated  by  arrows  on  the  abscissa;  as  shown  in  Fig.  14,  these  predicted 
FSP  locations  are  in  excellent  agreement  with  the  experimental  locations.  The 
experimental  Cw  values  become  significant  at  an  axial  station  about  one  duct  radius 
upstream  of  the  corresponding  FSP.  Indeed,  the  experimental  distribution  for  Ct  = 0.137 
indicates  that  pure  secondary-stream  fluid  exists  nowhere  in  the  flow  field,  even  at  x = 0. 

The  predicted  C^  distributions  in  Fig.  15  show  that  the  wall  concentration  remains 
zero  for  some  distance  downstream  of  the  FSP,  then  rises  rapidly  to  the  fully  mixed 
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value.  The  basic  integral  analysis  was  formulated  with  the  assumption  that  the  turbulent 
Schmidt  number  (Set)  is  unity.  Computations  for  the  Ct  = 0.345  flow  were  also  made 
with  Set  = 0.7  in  the  wake  regimes;  as  shown  by  the  dashed  line  in  Fig.  15,  this  increase 
in  the  eddy  diffusivity  causes  the  predicted  Cw  to  begin  to  rise  at  the  FSP.  But  even  with 
Set  = 0.7,  the  predicted  distribution  is  in  poor  agreement  with  the  experiment. 
Clearly,  the  axial  diffusion  of  jet  species  upstream  of  the  FSP  is  significant;  of  course, 
this  streamwise  turbulent  diffusion  is  neglected  in  the  basic  integral  analysis.  The  best 
that  can  be  said  of  the  predicted  C^  distributions  is  that  the  fully  mixed  value  is  reached 
at  approximately  the  correct  axial  station. 


Figure  15.  Wall  concentration  distribution  for  the  experiments  of 
Becker,  et  al. 

The  efficiency  of  a recirculation  zone  as  a stirring  device  is  strikingly  demonstrated 
by  the  results  shown  in  Fig.  15.  For  this  geometry,  reattachment  does  not  begin  until 
x/rw  is  greater  than  four;  therefore,  the  velocity  field  is  very  nonuniform  for  the  axial 
distances  shown  in  Fig.  15.  On  the  other  hand,  the  concentration  profiles  become 
uniform  (Cw  = Cc)  at  an  axial  station  only  about  one  duct  radius  downstream  of  the 
FSP.  Of  course,  this  observation  is  valid  only  for  small  values  of  rn/r^ ; for  large  values  of 
rn/rw,  such  as  investigated  by  Heskestad,  the  concentration  profiles  remain  nonuniform  at 
distances  far  downstream  in  the  third  regime. 
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3.4  EXPERIMENTS  OF  CHEDAILLE,  ET  AL. 

Chedaille,  Leuckel,  and  Chesters  (Ref.  20)  investigated  an  axisymmetric  recirculating 
flow  configuration  that  was  similar  to  the  one  investigated  by  Becker,  et  al.;  rn/rw  was 
0.031.  However,  Chedaille,  et  al.,  heated  the  jet  to  pve  Toji/Toai  = 1-51  and  made 
temperature  measurements  in  the  flow  field.  These  temperature  measurements  permit  tht 
distribution  of  jet  species  concentration  in  the  flow  field  to  be  deduced. 

Axial  distributions  of  wall  and  centerline  concentration  are  shown  in  Figs.  16a  and  b 
for  two  of  the  flows  investigated  by  Chedaille,  et  al.  Even  though  these  flows  involve 
density  gradients,  the  experimental  results  are  similar  to  those  obtained  by  Becker,  et  al. 


a.  u,i/Uji  = 0.0086,  Toj/To*  = 1.51 
Figure  16.  Wall  and  centerline  concentration  distributions  for  li.'* 
experiments  of  Chedaille.  et  al. 

The  Cw  distribution  begins  to  build  up  at  an  axial  station  well  upstream  of  the  FSP. 
However,  C^  and  Cc  do  not  approach  the  fully  mixed  values  as  rapidly  as  for 
constant-density  flow.  The  predicted  results  are  qualitatively  similar  to  those  obtained  for 
the  Becker,  et  al.,  geometry.  For  Set  = 1.  Qy  remains  zero  for  some  distance  downsteam 
of  the  FSP,  then  rapidly  rises  to  the  fuUy  mixed  value.  For  Set  = 0-7,  the  predicted  Cw 
distribution  begins  to  rise  just  downstream  of  the  FSP. 
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Axial  Distance,  x/r^ 

b.  Uai/Uji  = 0.0050,  Toj/Toa  = 1.51 
Figure  16.  Concluded. 

3.5  COMPARISON  WITH  OTHER  PREDICTIONS 

In  the  limiting  case  of  constant-density  flow  with  a small  value  of  rn/r^  (nearly  a 
point-source  jet),  the  present  integral  analysis  becomes  similar  to  that  of  Hill  (Ref.  11). 
The  extent  of  recirculation  predicted  by  the  present  analysis  for  the  configuration  of 
BarchUon  and  Curtet  is  shown  in  Fig.  17,  along  with  HUl's  prediction.  The  two 
predictions  for  the  location  of  the  FSP  are  nearly  the  same.  However,  Hill's  analysis, 
which  is  based  on  a local  model  for  the  eddy  viscosity,  greatly  overpredicts  the  distance 
to  the  RSP.  The  present  theory  would  yield  similar  results  if  a local,  rather  than  a frozen, 
eddy  viscosity  were  used  in  the  reattachment  region. 

The  finite-difference  elliptic  analysis  of  Gosman,  et  al.  (Ref.  8),  was  applied  to  the 
flow  configuration  of  Barchilon  and  Curtet  (Ct  = 0.305).  Lateral  variations  of  the  eddy 
viscosity  were  assumed  to  be  negli^ble,  and  the  axial  distribution  of  the  eddy  viscosity 
was  prescribed  to  be  the  same  as  that  obtained  in  the  basic  integral  analysis.  Therefore, 
differences  between  the  elliptic  and  integral  solutions  are  attributable  to  the  different 
governing  equations  and  numerical  procedures,  rather  than  to  different  turbulent 
transport  properties. 

The  axial  distribution  of  centerline  velocity  predicted  by  the  elliptic  analysis  (Fig. 
18a)  is  in  better  agreement  with  experiment  than  is  the  prediction  of  the  integral 
analysis.  But  the  prediction  of  the  integral  analysis  for  the  axial  distribution  of  Uw  is 
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Figure  17.  Comparison  of  present  theory  with  Hill's  results  on 
extent  of  recirculation. 


better  than  that  of  the  elliptic  analysis  (Fig.  1 8b).  The  elliptic  prediction  of  the  location 
of  the  FSP  is  particularly  poor.  This  result  must  be  attributed  to  the  numerical 
procedures  or  to  the  specification  of  the  vorticity  boundary  conditions,  because  the  flow 
upstream  of  the  FSP  is  nearly  parabolic  and  the  eddy  viscosity  model  used  is  realistic  in 
this  region.  (Experience  has  shown  that  the  location  of  the  recirculation  region  predicted 
by  the  elliptic  analysis  can  be  significantly  affected  simply  by  changing  either  the 
computational  mesh  or  the  method  of  specifying  the  vorticity  at  the  primary  nozzle  lip 
(Refs.  9 and  17.)  In  the  present  computations,  these  factors  were  selected  after  some  trial 
and  error  to  yield  a good  prediction  of  the  centerline  velocity  distribution.) 

In  the  computations,  the  Barchilon  and  Curtet  flow  was  assumed  to  have  tracer 
species  in  the  primary  fluid.  Predicted  axial  distributions  of  the  concentration  of  this 
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Axial  Distance,  x/r^ 


a.  Centerline  velocity  distribution 


b.  Velocity  near  wall 

Figure  18.  Comparison  of  integral  theory  and  elliptic  theory 
for  the  Barchilon  and  Curtet  geometry. 
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c.  Wall  and  centerline  concentration 
Figure  18.  Concluded. 

tracer  along  the  wall  and  along  the  centerline  are  shown  in  Fig.  18c.  The  elliptic  solution 
yields  the  physically  realistic  result  that  Cw  begins  to  increase  upstream  of  the  FSP. 
However,  with  the  eddy  viscosity  model  used,  the  computed  counterstream  diffusion  of 
species  along  the  wall  is  less  than  shown  in  the  experiments  of  Becker,  et  al.  (Note  that 
the  fully  mixed  concentration  predicted  by  the  elliptic  analysis  is  higher  than  the  correct 
value  of  0.250.  However,  in  this  particular  computation,  the  error  is  only  about  10 
percent.) 

One  of  the  major  objectives  of  making  the  calculations  with  the  Gosman  program 
was  to  evaluate  the  species  profile  shape  function  that  is  used  in  the  integral  analysis.  The 
relationship  in  the  shear  layer  between  the  normalized  concentration  and  the  normalized 
velocity  is  shown  in  Figs.  19a  through  c.  The  integral  and  elliptic  results  are  compared  at 
equivalent  values  of  u^,  rather  that  at  equivalent  axial  stations.  For  flows  without 
recirculation,  the  integral  analysis  predicts  that  the  concentration  profile  parameter  (K) 
remains  nearly  one  (Fig.  19a).  However,  downstream  of  the  FSP,  K increases  rapidly  as 
the  recirculation  velocities  become  larger  (Figs.  19b  and  c).  The  elliptic  analysis  also 
predicts  a large  distortion  of  the  concentration  profile  relative  to  the  velocity  profile.  So, 
even  through  the  concentration  profiles  predicted  by  the  two  analyses  are  different  in 
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c.  Uw/Uji  = -0.018 
Figure  19.  Concluded. 

detail,  one  can  conclude  that  the  large  distortion  of  the  concentration  profile  predicted 
by  the  integral  analysis  is  qualitatively  correct. 

3.6  DISCUSSION 

The  integral  analysis  shows  that,  for  constant-density  axisymmetric  flows,  the  Curtet 
number  is  a very  good  correlating  parameter  for  the  occurrence  of  recirculation,  even  for 
rn/r*  values  much  larger  than  are  usually  considered.  Predicted  limiting  values  of  the 
initial  velocity  ratio  (uai/uji)  that  just  eliminate  recirculation  are  shown  in  Fig.  20.  Of 
course,  the  limiting  velocity  ratio  depends  strongly  on  r„/rw.  But  when  the  results  are 
interpreted  in_  terms  of  Q,  the  limiting  value  of  the  Curtet  number  is  seen  to  change  by 
only  10  percent  as  rn/rw  is  increased  from  0.02  to  0.50  (Fig.  21). 

Based  on  the  comparisons  with  experiment  that  have  been  presented,  several 
conclusions  can  be  drawn  about  the  basic  integral  analysis  of  ducted  flows  with 
recirculation.  First,  the  overall  velocity  field,  including  the  location  of  the  front  and  rear 
stagnation  points,  is  fairly  well  predicted,  at  least  for  flows  with  small  density  gradients. 
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The  major  deficiency  in  the  velocity  predictions  is  the  unrealistic  cusp  in  the  wall 
velocity  distribution.  The  predicted  concentration  field  is  much  less  satisfactory  than  the 
predicted  velocity  field;  a major  deficiency  in  the  analysis  is  the  neglect  of  counterstream 
turbulent  diffusion  of  species  near  the  FSP.  Note  that,  for  the  flows  that  have  been 
considered,  the  density  is  uncoupled  or  is  only  weakly  coupled  to  the  concentration. 
Therefore,  errors  in  the  concentration  field  have  only  a small  influence  on  the  predicted 
velocity  field.  However,  when  the  density  is  strongly  coupled  with  the  concentration,  as 
it  is  in  the  hydrogen-air  experiments  of  Schulz  (Ref.  9)  and  Chriss  (Ref.  21),  the 
erroneously  predicted  concentration  field  can  be  expected  to  lead  to  significant  errors  in 
the  predicted  velocity  field. 


Figure  20.  Limiting  velocity  ratio  for  recirculation  in  a 

constant-density  flow  within  a cylindrical  duct. 
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Clearly,  the  major  deficiencies  of  the  basic  integral  analysis  are  the  neglect  of 
streamwise  diffusion  of  species  (which  is  particularly  important  near  the  FSP)  and  the 
neglect  of  streamwise  turbulent  transport  of  momentum  (which  is  particularly  important 
near  the  beginning  of  reattachment).  A deficiency  of  less  importance  is  the  neglect  of  the 
small  but  significant  axial  pressure  gradients  that  occur  between  the  FSP  and  the 
beginning  of  reattachment.  An  extended  integral  analysis,  in  which  these  deficiencies  are 
avoided,  is  formulated  in  Section  4.0. 


Figure  21.  Limiting  Curtet  number  for  recirculation  in  a 
constant-density  flow  within  a cylindrical  duct. 

4.0  FORMULATION  OF  EXTENDED  INTEGRAL  ANALYSIS 

All  of  the  assumptions  used  in  formulating  the  basic  analysis  are  retained  in  the 
extended  analysis,  with  two  exceptions.  First,  streamwise  turbulent  transport  is  not 
assumed  to  be  negligible  in  the  first  wake,  second  wake,  first  wall  and  third  regimes  (The 
analysis  of  the  first  and  second  regimes  is  unchanged).  Second,  it  is  not  assumed  that  the 
static  pressure  is  constant  in  the  first  wake  and  second  wake  regimes;  instead,  the  axial 
distribution  of  static  pressure  in  these  regimes  will  be  computed. 

4.1  BASIC  EQUATIONS 

The  basic  differential  equations  that  describe  the  flow  field  are: 

Continuity 

(pur®)  ( _<i  (pvr“)  =0  (1) 

d\  i9r 
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Axial  Momentum 


pur“  ^ + pvt“  ^ = — r“  ^ (rr“)  + (pr®)  (32) 

dx  dr  dx  dr  dx 


where  r is  the  turbulent  shear  stress  and  a is  the  turbulent  normal  stress  in  the  axial 
direction. 

Element  Species  Conservation 

pur®  ^ + pvT®  ^ (qr®)  + _i(j8r“)  (33) 

dx  dr  dr  dr 


where  q is  the  lateral  turbulent  species  flux  and  is  the  turbulent  species  flux  in  the 
axial  direction. 

As  in  Section  2.3,  Eqs.  (1),  (32).  and  (33)  are  integrated  to  obtain  five  integral 
equations; 
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Half-Radius  Species 

q r“  + p u (C  — C ) r®  E 

m “c  c ' c dx 

r 

+ (37) 

•*  5x 

r 

c 

Equations  (34)  through  (37)  are  the  same  as  Eqs.  (6)  through  (9),  except  for  the  last 
term  that  appears  on  the  right-hand  side  of  each  equation. 

The  static  pressure  is  no  longer  assumed  constant  in  the  first-wake  and  second-wake 
regimes.  Therefore,  an  additional  equation  must  be  introduced  to  permit  the  pressure  to 
be  calculated.  The  additional  equation  is  chosen  to  be  an  axial  momentum  equation  for 
the  flow  near  the  wall.  Because  t and  3u/9r  are  assumed  to  be  zero  in  this  region,  the 
wall  momentum  equation  is: 


r J^(puC)r®dr  - C J _^(pu)r®dr  = 

J dx  dx 

r r 


du 

p U YL  = 


(38) 


4.2  STREAMWISE  TRANSPORT  TERMS 


A gradient  model  of  the  streamwise  transport  processes  is  assumed.  Therefore, 


a - 


and 


/3  = pe^ 

dx 


(39) 


(40) 


where  the  eddy  viscosity  is  assumed  to  be  a function  of  only  x;  the  eddy  viscosity  model 
is  the  same  as  that  used  in  the  basic  integral  analysis  (Section  2.6). 


With  the  gradient  model  for  a and  |3,  the  streamwise  transport  terms  in  Eqs.  (34) 
through  (38)  contain  second  derivatives  of  u and  C with  respect  to  x.  These  second 
derivatives  cause  the  system  of  equations  to  be  elliptic  in  nature;  of  course,  the  solution 
techniques  used  for  solving  the  integrated  boundary  layer  equations  are  not  applicable  to 
this  elliptic  system.  However,  the  solution  techniques  described  in  Sections  2.7  and  2.8 
can  be  used  if  the  streamwise  transport  terms  are  computed  approximately. 

Leibnitz'  rule  is  applied  to  the  streamwise  transport  integral  terms  that  appear  on 
the  right-hand  sides  of  Eqs.  (34)  through  (37); 
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(41) 


(42) 


(43) 


(44) 


Now,  we  introduce  the  approximation  that,  at  any  x. 


-i-  /^  ) r“dr  = _L  )r“dr  - _1_ 
dxj  AxJ  Ax 


J'or^Ar 


(45) 


where  Ax  is  the  computation  step  size  and  the  subscript  denotes  a quantity  evaluated 
at  the  last  upstream  computation  station,  (x  - Ax).  With  Eq.  (45)  and  the  definitions  of  a 
and  |3,  Eqs.  (41)  through  (44)  become: 
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(48) 


(49) 


Similarly,  the  streamwise  transport  term  on  the  right-hand  side  of  Eq.  (38)  is 
approximated  by 


In  the  first-wake,  second-wake,  and  first-wall  regimes,  the  terms  in  Eqs.  (47)  and 
(49)  that  contain  drm/dx  are  nonlinear  in  the  derivatives  of  the  dependent  variables. 
Therefore,  in  these  regimes,  drm  /dx  is  approximated  by  (drm  /dx)g. 

4.3  METHOD  OF  SOLUTION 

With  the  approximate  streamwise  transport  terms,  Eqs.  (5),  (34),  (j:;.,  (36),  (37), 
and  (38)  can  be  transformed  into  a system  of  ordinary  differential  equations  that  is  linear 
in  the  derivatives  of  the  dependent  variables.  The  transformation  procedure  is  the  same  as 
used  for  the  basic  integral  analysis. 

First-wake  regime  - The  flow  is  completely  described  by  six  variables  (u^,  r,,  b,  C^, 
K,  p).  The  resulting  system  of  equations  is 

N|  N,  lla  - N,  + N.  ^ ^ N-  M f-  N,  if  = N-  (51) 

dx  ^clx  ^dx  ^dx  ^dx  °dx  ' 

where  N = F*  for  the  continuity  equation,  N = G*  for  the  overall  momentum  equation, 
N = H*  for  the  half-radius  momentum  equation,  N = I*  for  the  overall  species  equation, 
N = J*  for  the  half-radius  species  equation,  and  N = L*  for  the  wall  momentum 
equation.  The  equations  of  the  coefficients  for  all  the  regimes  are  presented  in  Appendix 

B. 
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Second-wake  regime  - The  six  dependent  variables  are  Uc,  b,  C*,  K,  p),  and  the 
resulting  system  of  equations  is 


N,  —sr 

* tlx 


<li  ..  It 

>0  — + N, 


iJC 


d\ 


(ix 


N.  + N. 


dx 


dx 


N, 


(52) 


Third  regime  - The  waU  momentum  equation  is  not  used  in  the  third-  or  first-wall 
regimes.  The  five  dependent  variables  are  (u*,  Uc,  p,  C^,  K)  and  the  system  of  equations 
is 


<)  U 


dC 


dx 


- dx 


(53) 


fi  X 


dx 


dx 


First-wall  regime  - The  five  dependent  variables  are  (u*,  ri,  p,  Cw,  K),  and  the 
system  of  equations  is 


N,  N,  ^ + N,  ^ + N.  ^ + N-  dK  = N 


dr 


dC 


(54) 


dx 


dx 


dx 


dx 


With  the  coefficients  evaluated  numerically,  Eqs.  (51)  through  (54)  are  solved  with 
the  same  techniques  that  are  used  in  the  basic  integral  analysis. 

4.4  INITIATION  OF  THE  WAKE  REGIMES 

In  order  to  account  for  the  counterstream  turbulent  diffusion  of  species  in  the 
vicinity  of  the  FSP,  the  first-wake  or  second-wake  regime  (whichever  is  appropriate  for 
the  particular  flow)  is  initiated  upstream  of  the  FSP.  Thus,  the  wake  regime  that  is  first 
encountered  in  a given  flow  is  initiated  with  a positive  value  of  Uw , rather  than  with  Uw 
= 0 as  in  the  basic  analysis. 

As  the  outer  flow  approaches  the  FSP,  the  velocity  decreases  toward  zero.  In  a 
region  just  upstream  of  the  FSP,  the  velocity  is  sufficiently  low  so  that  the  turbulent 
motion  in  the  recirculation  region  can  propagate  upstream.  Thus,  the  appropriate  wake 
regime  must  be  initiated  at  the  axial  location  where  the  approach  velocity  (ua ) is  equal  to 
an  average  turbulent  fluctuation  velocity  (u')  in  the  outer  flow. 

The  fluctuating  velocity  (u')  is  calculated  by  assuming  that  the  eddy  viscosity  in  the 
outer  flow  is  the  same  as  in  the  shear  layer  at  the  particular  axial  station.  But  the  eddy 
viscosity  (e)  is  the  product  of  a turbulent  velocity  scale  and  a turbulent  length  scale  (£): 

e = (55) 
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The  length  scale  in  the  outer  flow  is  assumed  to  be  proportional  to  the  width  of  the 
outer  flow,  that  is,  to  [r*  - (ri  + b)].  With  the  eddy  viscosity  defined  by  Eq.  (23),  the 
turbulent  velocity  scale  in  the  outer  flow  is  described  by 

(56) 

= ['w 

where  c is  a prescribed  constant,  which  should  be  about  one.  At  each  axial  station  in  the 
first  and  second  regime  computations,  Ua  is  compared  with  u'  obtained  from  Eq.  (56). 
The  appropriate  wake  regime  is  initiated  when  Ua  = u'. 

4.5  DISCUSSION  OF  THE  EXTENDED  ANALYSIS 

Although  streamwise  transport  terms  have  been  included  in  the  momentum  and 
species  conservation  equations,  the  analysis  must  be  considered  at  best  only 
"quasi-elliptic,"  since  it  was  developed  only  to  compute  approximately  the  locally  elliptic 
regions  in  flows  that  are  grossly  parabolic.  Because  of  the  assumptions  in  the  analysis,  the 
downstream  boundary  conditions  are  implied  to  be  only  weakly  coupled  to  the 
recirculation  region.  That  is,  the  duct  must  have  no  abrupt  changes  in  cross-sectional  area, 
and  must  be  long  enough  so  that  the  flow  downstream  of  the  RSP  approaches  the  fully 
mixed  condition.  In  essence,  the  downstream  boundary  conditions  are  assumed  to  be 
isolated  from  the  elliptic  region  by  a region  of  parabolic  flow. 

5.0  CONCLUDING  REMARKS 

The  basic  integral  analysis^  which  involves  the  solution  of  the  integral  form  of  the 
boundary-layer  equations,  is  computationally  efficient  and  yields  results  that  are 
sufficiently  accurate  for  many  purposes.  Indeed,  the  predictions  for  the  location  and 
extent  of  recirculation  in  ducted  flows  are  very  good,  at  least  for  flows  with  moderate 
density  gradients.  But  the  streamwise  turbulent  transport  terms,  which  are  neglected  in 
the  basic  analysis,  are  important  in  establishing  the  concentration  Held  near  the  front 
stagnation  point  and  in  establishing  the  velocity  field  near  the  onset  of  reattachment. 

In  the  extended  integral  analysis,  which  has  been  programmed  for  numerical 
solution,  the  streamwise  transport  terms  have  been  included  in  an  approximate  manner. 
The  validity  of  this  quasi-elliptic  analysis  must  be  established  by  comparison  with  the 
available  experimental  results  on  constant-density  and  variable-density  ducted  flows  with 
recirculation. 

The  eddy  viscosity  model  used  in  this  study  yields  satisfactory  results  for  the  flows 
considered.  But  for  complex  variable-density  flows,  such  as  occur  in  combustion  systems, 
the  simple  eddy  viscosity  model  can  be  expected  to  be  unsatisfactory.  A more 
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sophisticated  turbulence  model,  based  on  solution  of  the  turbulent  kinetic  energy  (TKE) 
equation,  will  be  required  for  complex  flows.  For  flows  without  recirculation,  the  present 
integral  analysis  has  already  been  extended  to  include  solution  of  an  integrated  TKE 
eqution  for  the  shear  layer  (Ref.  22). 
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APPENDIX  A 

COEFFICIENTS  FOR  THE  BASIC  INTEGRAL  ANALYSIS 

The  following  parameters  are  defined : 

y = (r  - r;)/b 

r^dr  = ab^ydy  + r®  bdy  ; a = 0 or  I 
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1 


z 


k 


-d—  (pit) 
i9z. 
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2z 


k 


-d—  (pu^) 
dz^ 


(puC) 


" 5fr‘"' 


S 


5z, 


(C) 


where  is  a generalized  variable. 

The  first  and  second  regime  coefficients  (F,  G,  H,  and  I)  are: 


F,  ^ ab2 


J Su,  yjy  + '?b  J s„,^  dy  + ^ (p,"|) 

o o 

r r«+*  - (r.  + 1))“+’  "I 

" [—  ;rn  J ^ 
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= ab2 


1 I 


dy 


f's  = 


1 1 
Jsj^^ydy  - r“b  Jsj^^dy  - r« 
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dx  dx 
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j^2z^  ydy  + r<fb  dy 


Gs  = 


-ah^ 


jhz^y^y  - js^.^dy  - p,u2  r«  r« 


r“  r 
w w 


llj  = ab' 


1/2  1/2 
J Sa.^ydy  * J 


S2j:  ^y  - "m  «'' 


1/2 

j >''*>■ 


Hg  = ab 


//  a-l  ,a+l\  _ „ , -,  /ja-*-l_  a+ 

S,  dv  -(-J E — ) f-i  (PiU?)  - ± (p.u,)1-(_^s ^ 

’='1  ■ \ a+1  / [dp  J J dp  ^ ' J \ “+' 

O 

1/2  1/2  1/2  1/2 
^ J ^2,^  ydy  + r?b  J Sj^^dy-  u^ab^y*  Sj^^ydy  - u^r?bj  Sj^^dy 


52 


AEDC-TR-77-115 


1/2  1/2  1/2  1/2 
\]^  = ab^  J ■‘’2-.,,  f ~ ^ f 


U,  = ab^ 


1:2  1/2 

J j 
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1 / 2 1/2 


II- 

o 


1/2  1/2  1/2 
■ ah^  J ^27^y^y  ~ j ^27..y'^y  + “m  J ^iz. 


1/2 


I,  = 


('"’.7^-)  IF  ‘'’••“i’  * '*'’' 


dv 


1 1 
h = p3z2>'‘l>-  + I’  '*>■ 


*3  “ 


I.  = ab' 


1 1 
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^3z  ^y 


I,  = -ab' 


y^3z.  ydy  - r^bj'sg^^dy  -r  PcU^Cj,r^“ 


In  the  first  regime  (where  Pc  = Pj,  Uc  = uj  and  Q = 1),  Zi  = p,  Z2  = rj,  zj  = b,  Z4  = K,  and  Z5 
= X.  In  the  second  regime  (where  ri  = rc),  zj  = p,  Z2  = Uc , Z3  = b,  Z4  = k,  and  Z5  = x. 
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The  third,  first  wake,  second  wake  and  first  wall  regime  coefficients  (F,  G,  H,  1,  and 
J)  are: 


F,  . ab2  * [' * ' ^ ] 
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1/2 


1/2  1/2 
y*  dy  + Cn,  ab2  y S^^^ydy 

1/2 

+ C„,r?'b  J Si,,^dy  + Pc«c 


_ p ) r ^ — ■!“  Q r 


In  the  third  regime  (where  rj=  Tc,  b = r^  - Tc,  and  6 = 1),  zj  = , Z2  = Ug,  Z3  = p,  24  = 

Cw , Z5  = K,  and  Zg  = x. 
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In  the  first-wake  regime,  zi  = , Z2  = rj,  Z3  = b,  Z4  = C* , Z5  = K,  and  zs  = x (the 

pressure  is  constant  in  this  regime;  therefore  8 = 0). 

In  the  second-wake  regime  (where  = rg),  zi  = u*,  Z2  = Uc,  Z3  = b,  Z4  = C^,  Z5  = K, 
and  Z5  = X (the  pressure  is  constant  in  this  regime;  therefore  S = 0). 

In  the  first-wall  regime  (where  b = - rj  and  S = 1 ),  zj  = , Z2  = Tj,  Z3  = p,  Z4  = , 

Z5  = K,  and  zg  = x. 
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APPENDIX  B 

COEFFICIENTS  FOR  THE  EXTENDED  INTEGRAL  ANALYSIS 


The  parameters  defined  in  Appendix  A are  also  used  in  this  appendix.  Whenever 
possible,  the  coefficients  are  written  as  the  basic  analysis  coefficients  with  additional 
terms  as  required. 

In  the  first-wake  and  second-wake  regimes,  the  coefficients  (F*,  G*,  H*,  I*,  J*,  L*) 

are: 
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dv  - e p^r“  ^ :!£c  + _ rC 
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= h - ih  I ^57.,  y^y  - 


1/2 

-i.  r<fb  f S. 


dv  — < p r“  E.  S 

r.„  ' c"c  c ^ 3 

2 d X d 11 


- ep  r®  (^'*'^"1^  " m 

\ dx  /g  d,. 


dr  \ do 


■I3*  = ‘*3  - ^ ^5z  y^y  - I ‘‘^sz  ‘ly 


1/2  1/2 

J,.  . J4  - ^ ah^J  S,,^yi,  - ^ -?by  S.,^dy  - ,p„r« 


Js*  - -Is 


1/2  1, 


Ax  ■ I ="5 


- f p r®  ^l£  + fP  '■'*  (dr  /dx)n 

dx  dK  m ^ m'  'g 


1/2  1/2  1/2 

V = J y^y  + j J y^y 

00  o 

1/2  1/2  1/2 

- C„r?b  f S,  dy  - _«_  ab2  f S,  ydy  - _L  r?b  / S. 

/ 1^6  Ax  / ^ J « 

00  O 

r,a+l  _ ,a+n  _ r®+n 

+ ij s — i (PiU;)  - C Li £_J  J-  (pjUj) 

(a  +1)  dp  J ■"  (a+  1)  dp  J J 


- e„p„T°_  ^ + fp„r®  (dr  /dx)n 

‘=^'=  '=  dx  dp  * dp 


T * - 1 - 1 
''7  - h - XZ 


* m *■  f ^ *• 

- h f ^ * fx  ^ / '5.,  <*>' 

J o o 

-Ig 


Li*  = p u — _£.  p 
1 r\H  V!  A„  r'\ 


Lg*  = 0 
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Lg*  = 0 


L4*  = 0 


Lg*  = 0 


Lg-  = 1.0 


L-*  = -_L  fp  ^ 

Ax  dxj, 


In  the  first-wake  regime  (where  Q = 1 and  ec  = 0),  Zi  = , Z2  = n,  Z3  = b,  Z4  = , 

Z5  = K,  zg  = p,  and  Z7  = x.  In  addition, 

dn/dx  = 0 

dc  ^ dc  = dc  = dc  ^ n 

dC  dC  dC 

= m — m_  _ Q 

i9u  (?u  (9p 

w c * 

In  the  second-wake  regime  (where  ri  = rc  and  €c  = e),  zj  = u^.,  Z2  = Uc,  Z3  = b,  Z4  = 
Cw , Z5  = K,  zg  = p,  and  z^  = x. 

The  coefficients  (F*,  G*,  H*,  I*,  J*)  in  the  third-  and  first-wall  regime  are: 

F * - F 

1 ] - 1 j 

F * - F 
*^2  - *^2 

F * - F 
.3  “ 

F,*  = F, 

^5*  = Fg 
Ffi*  = Fg 


G,*  = G,  - 


1 1 

ah^  I S.  ydy  - r?b  j S.  dy 
Ax  J '‘zi  ^ ^ ' Ax  J ‘‘^1  ^ 


n <*>■ 
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= G„ 


1 1 

- t.  j 


?b  _f_  Is.  dv  — f p r“ 

1 A..  I 4Z2  - ^cPc  c ^ 


Ax 


s- 


Gg*  = Gs  - T*-  ab2  I $4^  ydy  - ^ r«b  | dy 


A} 


,.Js 


4z, 


G4*  = C4 


G5*  = G5 


Ax 


1 

1 

— 

r __ 

•w 

Js,,^ydy  . ^ ,?b 

j®-.  - h 

O'  dr 

J 

0 

J 

0 

:- 

Hj*  = H,  - ab^  ^ 


Ax 


1/2  y2 

j S4X,  ydy  - ± j S4,.^ 


dv  + 1.  ( p SI 

y + 2 dx 


1/2 


Hj*  = Hg  - ab- 


1/2 

I S^,  ydv  - rfb  _f_  f S.  dy  - 1 p r«  ^ -e  p r«  *!!£ 

Ax  J ^==2  ^ ' ‘ Ax  J *'‘2  ' 2 dx  ^ dx 


Hg*  = Hg 


1/2  1/2 


1 - „ dr  da 

dy  - f p_r„  — !2.  — £ 
J f'm  m 


H4*  = H4 


H5*  = Hg 


1/2 


1/2 


He* 


He  * i;-'’'  j y^y  "*'?■>/ 


S4/,  dy  - 1 


Ax 


J 


dr 


J S. 


Ij*  = Ij  - ^ ab' 


Ax 


1 1 

ydy  - -^-  r^b  Is,  dy  — f p r“  l!!® 

5*1  ^ ^ A*  J '1  d*  da^ 
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I * _ I 
Ig  - I2 


1 1 

ab^  I Sc  ydy  — r?b  | Sc  dy  — f p r“  S.  «= 

Ax  J ^'2  ^ J 2 ^ dx 


I,*  = I,  - -f- 

3 3 Ax 


ab^  jS-  ydy  - _f_  r?b  | ! 

J Ax  ■ J 


- ^ r“b  I Sc„  dy  - e„p„r“  ^ 

dx  (5p 


'4*  = I4 


1 1 

- -L  ab^  I Sc_  vdy  - _i_  r?b  |S-  dy 

Ax  J ' Ax  ' J ”1 


'5“ 


fp  r“  !!1«  - 6 p r®  ^ ^ 

"■  dx  •=  dx  dC. 


a 

^ f D r”  £ 

*=  dx  dK 


1 1 r r*  " 

-i-  ab^  I Sr,,  ydy  + -L.  r?b  | S-  dy  - .1=  I j8r“dr 

Ax  J ^ Ax  ' J ^^6  Ax  J 


'6*  = >6  + 


1 1 


, „ dr  dC  „ dr  dC 

dx  da  '=  dx  da 


J2*  - J2  - 


1/2 


* yd, 


1/2 

/• 


-L.  r?b  I S-  dv 
^2  ' 


a dr  (9C 

in_  y 

” m ni  J „ ^ ^ 


.a  K 

dx  du 


J3’  = 13 


1/2 

-i'7 


^sz.  y^y  - ^ / s-„  dy 


Ax 


1/2 

/ 


dr  dC 
.a  ■n  m 
m 


^ ^Pm  *'m  37  ^TT  ~ ^cP"  ’’ 


dx  dp 


dr  dC_ 
a c c 

~dr  dp 


J4  = J4  - 


1/2  1/2 

i;  “>■'  y S5.,  xJy  -f'^j 

0 0 


dy 


dr„  dC. 


+ <Pm‘‘ 


dx  dC.. 


^cPc*^c 


dC„ 


dx  dC. 
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J_  r“b 

Ax  ' 


1,'2 

/ 


dv 


fPm 


«lr 

m 

(lx 


dCm 

IT) 


dr 

^r  — 
' ilx 


dC^ 

dK 


- h 


1.-2  1:2 

* f.  “*■'  / ■‘‘w..  >■'*>■  * j 


S^,  dy  - ^ 

j Ax 


r 

/ 


L^i 


j9  r^dr 


In  the  third  regime  (where  r,  = rc  and  b = - rc),  zi  = , Z2  = Ug,  Z3  = p,  Z4  = Cw , Z5 

= K,  and  = x.  In  addition,  6^  = 6 and  3u/3p  = 3Uc/3p  = 0. 


In  the  first-wall  regime  (where  b = - rj,  Q = 1,  and  = 0),  zj  = u^,  Z2  = n,  Z3  = 

p,  7.4  = Cw , Z5  = K,  and  zg  = x.  In  addition. 


<9f: 

= !^  = 

0 

(9« 

w 

r?p 

r}(: 

r7(: 

3(: 

m 

.=  n)  _ 

in 

rtw 

t' 

(5r 

1 

3p 
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NOMENCLATURE 

b Shear  layer  thickness 

C Element  fraction  of  primary  stream  gas 

Ct  Curtet  number 

c Constant  in  Eq.  (56) 

Cfw  Wall  skin  friction  coefficient 

Cp  Specific  heat  at  constant  pressure 

F,G,H,I,J  Coefficients  in  system  of  ordinary  differential  equations  for  basic  analysis 

Coefficients  in  system  of  ordinary  differential  equations  for  extended 
I*,J*,L*  analysis 

FSP  Front  stagnation  point 

H Stagnation  enthalpy,  including  chemical  heats  of  formation 

k Eddy  viscosity  coefficient 

ko  Eddy  viscosity  coefficient  for  constant-density  flow 

K Parameter  in  species  profile  equation 

£ Turbulent  length  scale 

Mm  Mach  number  at  half-radius  control  surface 

N Generalized  coefficient  in  system  of  ordinary  differential  equations 

p Static  pressure 

q Lateral  turbulent  species  flux 

r Lateral  coordinate 

Tc  Lateral  distance  from  centerline  to  centerbody  surface 

Ti  Lateral  distance  from  centerline  to  inner  edge  of  shear  layer 

rm  Lateral  distance  from  centerline  to  half-radius  control  surface 
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rn  Primary  nozzle  exit  radius 

rw  Lateral  distance  from  centerline  to  duct  wall 

R Gas  constant 

RSP  Rear  stagnation  point 

S Generalized  integrand  (Appendixes) 

Set  Turbulent  Schmidt  number 

To  Total  temperature 

u Axial  velocity  component 

Uk  Characteristic  velocity  defined  by  Eq.  (31a) 

u'  Turbulent  velocity  scale 

V Lateral  velocity  component 

X Axial  coordinate 

Ax  Integration  step  size 

y Dimensionless  mixing  zone  coordinate  (Appendixes) 

ziZ2...Zk  Generalized  dependent  variable  (Appendixes) 

a Geometric  parameter  (one  for  axisymmetric  flow,  zero  for  planar  flow) 

|3  Turbulent  species  flux  in  axial  direction 

5 Numerical  factor  used  in  the  appendixes  (5  = 0 where  pressure  is  constant,  5 
= 1 when  pressure  is  variable) 

6 Turbulent  eddy  viscosity 

p Density 

a Turbulent  normal  stress  in  axial  direction 

T Turbulent  shear  stress 
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SUBSCRIPTS 

1 Conditions  in  plane  of  primary  nozzle  exit 

a Conditions  in  isentropic  secondary  flow 

c Conditions  at  centerline  or  at  centerbody  surface 

j Conditions  in  isentropic  primary  flow 

£ Indicates  a quantity  evaluated  at  last  upstream  computation  station 

m Conditions  at  half-radius  control  surface 

w Conditions  at  duct  wall 
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